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- - - Abstract. Some Bianchi type I viscous fluid cosmological models with a variable 

' cosmological constant are investigated in which the expansion is considered only in 

. two direction i.e. one of the Hubble parameter {H\ = ^) is zero. The viscosity 

' coefiicient of bulk viscous fluid is assumed to be a power function of mass density 

whereas the coefficient of shear viscosity is considered as constant in flrst case 
whereas in other case it is taken as proportional to scale of expansion in the model. 
The cosmological constant A is found to be positive and is a decreasing function of 
Q,^ , time which is supported by results from recent supernovae la observations. Some 

■ physical and geometric properties of the models are also discussed. 

Keywords: Cosmology; Bianchi type I universe; viscous fluid models; variable 
^ . cosmological constant. 
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o 

T^lj- ' 1. Introduction 

o ■ 

O . Models with a dynamic cosmological term A(t) are becoming popu- 

lar as they solve the cosmological constant problem in a natural way. 



^J^' There are significant observational evidence for the detection of Ein- 

stein's cosmological constant, A or a component of material content 
of the universe that varies slowly with time and space to act like A. 
^ . Recent cosmological observations by High -Z Supernova Team and Su- 

H \ pernova Cosmological Project (Garnavich et al, 1998; Perlmutter et 

al, 1997, 1998, 1999; Riess et al, 1998; Schmidt et al, 1998) suggest 
the existence of a positive cosmological constant A with magnitude 
K{Gh/(?) « 10~^^^. These observations on magnitude and red-shift of 
type la supernova suggest that our universe may be an accelerating 
function of the cosmological density in the form of the cosmological A- 
term. Earlier researches on this topic are contained in Zeldovich (1968), 
Weinberg (1972), Dolgov (1983, 1990), Bertolami (1986), Ratra and 
Peebles (1988), Carroll, Press and Turner (1992). Some of the recent 
discussions on the cosmological constant "problem" and consequences 
on cosmology with a time-varying cosmological constant have been 
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discussed by Dolgov (1993,1997), Tsagas and Maartens (2000), Sahni 
and Starobinsky (2000), Peebles (2003), Padmanabhan (2003), Vish- 
wakarma (1999, 2000, 2001, 2002), and Pradhan et al. (2001, 2002, 
2003, 2004). This motivates us to study the cosmological models in 
which A varies with time. 

The distribution of matter can be satisfactorily described by a per- 
fect fluid due to the large scale distribution of galaxies in our universe. 
However, observed physical phenomena such as the large entropy per 
baryon and the remarkable degree of isotropy of the cosmic microwave 
background radiation, suggest analysis of dissipative effects in cosmol- 
ogy. Furthermore, there are several processes which are expected to give 
rise to viscous effects. These are the decoupling of neutrinos during the 
radiation era and the decoupling of radiation and matter during the 
recombination era. Bulk viscosity is associated with the GUT phase 
transition and string creation. Misner (1967, 1968) has studied the 
effect of viscosity on the evolution of cosmological models. The role 
of viscosity in cosmology has been investigated by Weinberg (1971). 
Nightingale (1973), Heller and Klimek (1975) have obtained a viscous 
universes without initial singularity. The model studied by Murphy 
(1973) possessed an interesting feature in which big bang type of singu- 
larity of infinite space-time curvature does not occur to be a finite past. 
However, the relationship assumed by Murphy between the viscosity 
coefficient and the matter density is not acceptable at large density. 
Roy and Prakash (1977) have investigated a plane symmetric cosmo- 
logical models representing a viscous fluid with free gravitational field 
of non-degenerate Petrov type 1 in which coefficient of shear viscosity 
is constant. Bali and Jain (1987, 1988) have obtained some expanding 
and shearing viscous fluid cosmological models in which coefficient of 
shear viscosity is proportional to rate of expansion in the model and 
the free gravitational field is Petrov type D and non-degenerate. The 
effect of bulk viscosity on the cosmological evolution has been inves- 
tigated by a number of authors in the framework of general theory 
of relativity (Padmanabhan and Chitre, 1987; Johri and Sudarshan, 
1988; Maartens, 1995; Zimdahl, 1996; Pradhan et al, 1997; Kalyani 
and Singh, 1997; Singh et al, 1998; Pradhan et al, 2001, 2002). This 
motivates to study cosmological bulk viscous fluid model. 

Recently Bali and Jain (1997) have investigated Bianchi type I vis- 
cous fluid cosmological models. Motivated by the situations discussed 
above, in this paper, we shall focus on the problem with varying cos- 
mological constant in presence of bulk and shear viscous fluid in an 
expanding universe. We do this by extending the work of Bali and Jain 
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(1997) by including varying cosmological constant and the coefficient 
of bulk viscosity as function of time. This paper is organized as follows: 
The metric and the field equations are presented in section 2. In section 
3, we deal with the solution of the field equations in presence of viscous 
fluid. In section 3.1, we consider the coefficient of shear viscosity (r/) as 
constant whereas in section 3.2, rj is taken as proportional to scale of 
expansion in the model. Section 4 includes the solution in absence of 
shear viscosity. In section 5, we have given the concluding remarks. 



2. The Metric and Field Equations 

We consider the Bianchi type I metric in the form 

where B and C are functions of t alone. This metric depicts the case 
when one of the Hubble parameters (here i^i = ^) is zero, i. e. the 
expansion is only in two directions. The kinematic parameters are then 
related as 6 = — 3cr|. This condition leads to the metric (1). 
The Einstein's field equations (in gravitational units c = 1, G = 1) 
read as 

E(-^Rgi + Agi = -87rT^, (2) 

where is the Ricci tensor; R = g^^ Rij is the Ricci scalar; and is 
the stress energy-tensor in the presence of bulk stress given by Landau 
and Lifshitz (1963) 

Ti = (p + p)viV^ + pg{ - r] (y{. + -i + v^v^Vi-i + ViV^v^.^ 

0{g\ + viv^). (3) 

Here p, p, rj and ^ are the energy density, isotropic pressure, coefficient 
of shear viscosity and bulk viscous coefficient respectively and vi is the 
flow vector satisfying the relations 

gijv'iP = -1. (4) 

The semicolon (; ) indicates covariant differentiation. We choose the 
coordinates to be comoving, so that 

= Q = v'^ = v^,v'^ = 1. (5) 
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The Einstein's field equations (2) for the fine element (1) has been set 
up as 

-B44 , C44 B4C4 



-87r 



P 



+ 



C 



+ 



BC 



+ A, 



-Stt 



-Stt 



■ . B4 



B4C4 
BC 



Bu 



+ A, 



A, 



(6) 
(7) 
(8) 
(9) 



where the suflfix 4 at the symbols A and B denotes ordinary differ- 
entiation with respect to t and 6 is the scalar of expansion given 

by 

e = v'.i. (10) 



3. Solution of the Field Equations 

We have revisited the solutions obtained by Bali and Jain (1997). Equa- 
tions (6) - (9) are four independent equations in seven unknowns A, 
B, p, p, r), ^ and A. For the complete determinacy of the system, we 
need three extra conditions. The research on exact solutions is based 
on some physically reasonable restrictions used to simplify the Einstein 
equations. 

From Eqs. (6) - (8), we obtain 

-B44 B4C4 i« -^4 

^ + ^ = -16^r?-, (11) 

and 

-B44 C44 f B4 C4\ 

Equation (12) on integration leads to 

C2(^|)^ = Le-i6-'?*, (13) 

where L is a constant of integration. Setting BC = p and ^ = v in 
Eqs. (11), (12) and (13) lead to 

/X44 + 167r?7//4 = 0, (14) 

and 

= £^-iGnvt_ (15) 



purnimalsty.tex; 7/02/2008; 3:00; p. 4 



SOME BIANCHI TYPE I VISCOUS FLUID COSMOLOGICAL MODELS .... 

Equation (14) on integration leads to 

where M is an integrating constant. 
Here we consider two cases: 



3.1. CASE I : LET rj = CONSTANT = a (SAY) 
In this case Eq. (16) on integration gives 

M = iV - -|^e-i«-*, 
loTra 

where is a constant of integration. 
Equation (15) on integration leads to 

u = a{16TraN - Me~^^™*)^, 

where a is an integrating constant. 
Now we set 

-Wvat 



e " = cos 2Vl67rar, 



(167ra) M 



and 

N ^ 



IGvra 

Using the above Eqs. (19) - (21) in (17) and (18), we obtain 

„ , , /sin \/l67rar\ ^ 

and 

( — 
sinVl67rar\ ^ 

Prom above Eqs. (22) and (23), we obtain 

„9 ^ 14-^ / sin -v/lGvrar \ ^ 

B^=,.= (2Mr^ 



_ 2^ 



\/167ra y 
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Hence the metric (1) reduces to the form 



tan 2\/167rar 



i+— 



L I sm 



2-1- — 



7ra 



+(2M)^-M 



1 _ A / sin VlSvrar 



dz^ 



Tra 



(26) 



After suitable transformation of coordinates metric (26) reduces to the 
form 



ds^ = -4 



/tan2/cr 
I k 



2 /sinfer\2+M I / sin kT 



dT^+dX^+( 



dY^+ 



where k = \/l67ra. 

The pressure and density for the model (27) are given by 



dz^ 
(27) 



Sirp 



1 



4M2 VsinfcT 



M(M-L) sin^ fcT cos 2kT- 



(L + M) 



{{L + M) cos 2kT - 2M} + 167rM2 (^^ - cos 2A;T 
8^P=7T7^(l7rT77;) ^ ^ ^ cos^ 2fcr 



4 

sin kT\'^' 



cos2kTx 



4M2 Vsinitr 



k 



+ A 



-A 

(28) 

(29) 



For the specification of ^, we assume that the fluid obeys an equation 
of state of the form 

P = IP: (30) 

where 7(0 < 7 < 1) is constant. Thus, given ^{t) we can solve for the 
cosmological parameters. In most of the investigation involving bulk 
viscosity is assumed to be a simple power function of the energy density 
(Pavon, 1991; Maartens, 1995; Zimdahl, 1996) 



e(t)=eop", 



(31) 



where ^0 and n are constants. If n = 1, Equation (23) may correspond 
to a radiative fluid (Weinberg, 1972). However, more realistic models 
(Santos, 1985) are based on n lying in the regime < n < \. 
On using (31) in (28), we obtain 



87rp 



1 



k 



4M2 VsinfcT 



M{M-L) sin^ kT cos 2kT- 



{L + M) 



cos2kTx 
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sin kT 



{(L+M)cos2A;r-2M} + 167rM2 ^^op" - ^aj cos2kT ^ 



-A 
(32) 



3.1.1. Model I: Solution for ^ = 

When n = 0, Equation (31) reduces to ^ = = constant. Hence in 
this case Equation (32), with the use of (29) and (30), leads to 



647rM2(l + 7)p 



k 



sin kT 



cos 2kT 



M^-L^ + 2{L^ + M^) sin^ kT+ 



„^ ,.2 2a\ /sinA;r\ 

32rf ({„--) (— ) 



2-1 



EHminating p{t) between Eqs. (29) and (33), we have 



2(1 + 7)A = cos2A;r 



(1-7)(M2-L2)+ 



2{(3M^ + L^) + (M^ - L^)j} sin^ kT + 64ttM^ 6 " ^ 



(33) 



2a\ /sinfcry 



— 



(34) 



3.1.2. Model I: Solution for ^ = S,oP 

When n = 1, Equation (31) reduces to ^ = ^Qp. Hence in this case 
Equation (32), with the use of (29) and (30), leads to 



k'^ cos 2kT 



M2[2(l + 7) sin^ kT - fc^^o cos 2kT] sin^ kT 



MLJ^ + (l2 + m2) sin2 kT-^ cos 2kT 
Eliminating p{t) between (29) and (35), we have 



k 



sin kT 



A = - 



(M2 - l2) cos^ 2kT ( k 



16M2 



sinfcT 



+ 



2M2[2(1 + 7) sin^ kT - k'^^o cos 2kT] sin^ kT 



+ (L2 + M2)sin2 kT 



■ cos 2kT 



k 



sin kT 



(35) 



(36) 



Our study for these models shows constant value of cosmological con- 
stant (A) for large values of time and do not decrease with time (this 
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means that the universe is not expanding or may be steady state con- 
dition). In this case, detailed study shows that the scalar of expansion 
9 does not increase with time. Our study is inconsistent with work 
done by Bali and Jain (1997). Bali and Jain claim that the universe 
is expanding which does not match with our result. Also the claim of 
minimum and maximum expansion rate in 6 is reflection of periodicity 
of trigonometrical functions involved there. We are trying to find fea- 
sible interpretation and situations relevant to this case. Further study 
is in progress. 

Some Physical Aspects of the Models: 

With regard to the kinematic properties of the velocity vector in 
the metric (27), a straight forward calculation leads to the following 
expressions for the scalar of expansion {0) and for the shear (cr) of the 
fluid. 

The rotation ui is identically zero. The non-vanishing components of 
conformal curvature tensor are 

The models represent shearing, non-rotating and Petrov type I non- 
degenerate in general, in which the flow is geodetic. However, if L = 
then spacctimc reduces to Petrov type D. The model starts expanding 
at T > but the initial expansion is slow. When T is closer to 7r/2fc, 
it has stiff rise in the expansion and then decreases. This shows the 
case of oscillation. It is observed that the expansion is minimum at 
T = or T = vr/A;. The large values of 9 near T = ■K/2k is reflection of 
trigonometric property. But expansion remains finite. The expression 
for I and ^ are found to be 
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^2 327rM2 cos2 2kT 



(M^ - L^) cos^ 2kT + IGM^A 



sin kT'" 



(43) 

Prom Equation (42), it can be seen that shear is proportional to scalar 
of expansion 9 in the model. Since limT->cx3 f 7^ 0, hence the models do 
not approach isotropy for large values of T. 

The rate of expansion Hi (Hubble parameters) in the direction of X, 
Y, Z are given by 

Hi=0 (44) 



3.2. CASE II: LET r] = (39, WHERE /3 IS AN ARBITRARY 
CONSTANT. 

In this case Equation (14) reduces to 

fifi44 + I67rpl4 = 0, (47) 

which on integration leads to 

//= [(l + 167r/3)(A;it + fc2)]™^, (48) 

where ki and ^2 are constants of integration. 
Equation (12) reduces to 

m(^) =-{l + 16nP)^, (49) 
\ v / 4 fJ, 

which on integration leads to 

zy = A;3(A;it + A;2)'=i(i+i«'^«, (50) 

where k-^ and /C4 are constants of integration. 
Prom Equations (48) and (50), we obtain 

B^ = ^u = k^kzikit + A;2)*^6(fei+fc4)^ (51) 
C2 = ^ = ^{k^t + ^2)^=6(^1-^4)^ (52) 

where ^ 

k5 = {l + 167r/3) (1+16-/3) , (53) 



pumimalsty.tex; 7/02/2008; 3:00; p. 9 



10 



A. PRADHAN AND P. PANDEY 
1 



A;i(l + 167r/3)' 



(54) 



Hence the metric (1) reduces to 



(55) 

After suitable transformation of coordinates, the metric (55) takes the 
form 



ds 



The pressure and density for the model (56) are given by 

kjk6{ki + k4) 



S-rrp = 

SirkTke 
+ J; 



4^2 



-{kike + k^ke - 2) 



ST 



(/Ci - /C4) + e 



k2 



l,2h.2(u2 _ u2\ 



4r2 



On using (31) in (57), we obtain 

kjke{ki + k4) 



Sirp = — - 



4^2 



{kike + k^ke - 2) 



+ i 



ST 



(A^i-M+CoP^ 



2fc2y2fcife6 



A 



(56) 



(57) 
(58) 



(59) 



3.2.1. Mode/ /.• Solution for ^ = Co 

When n = 0, Equation (59), with the use of (58) and (30), leads to 



87r(l + 7)p 



klkQ{ki + k4){l - k4k6) 



+ ^ 



3T 



2r2 



{kl - ^4) + Co 



A;2 



Eliminating p(i) between Equations (58) and (60), we have 
kike{ki + k4 



(1 + 7)A 



_^87rA;f fee 



3r 



2 - {fci + fc4 + {kl - k4)j}k6 
k2 



{kl - k4) + Co 



2fe2T2fcife6 



(60) 



(61) 
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3.2.2. Model I: Solution for ^ = ^oP 

When n = 1, Equation (59), with the use of (58) and (30), leads to 



87r 



1 + 7 



6T2 



3(A:i + A:4)(l-A:4/c6) 



+167r/3A;i(fci - 3^4) 



Eliminating p{t) between (58) and (62), we have 



1 + 7- 



A = 



klkeih + ki) 

4r2 



2-{ki + ki)kfi 



(62) 



-(fci - fc4)A;6{7 



+ 



STrPkfk^iki - 3/^4) 
3T2 



(63) 



Prom Equations (61) and (63), we observe that the cosmological con- 
stant in both models is a decreasing function of time and it approaches 
a small value as time progresses (i. e. the present epoch), which explains 
the small value of A at present. 
Some Physical Aspects of the Models: 

With regard to the kinematic properties of the velocity vector in 
the metric (56), a straight forward calculation leads to the following 
expressions for the scalar of expansion (9) and for the shear (a) of the 
fluid. 



a 



T 

k-^ kQ {k 



2r2 



(64) 



(65) 



The rotation lo is identically zero. The non-vanishing components of 
conformal curvature tensor are 



'-^12 — 



^13 



i2r2 

12T2 

0]_4 



\ki — 3A;4 + A;4A;6(3A;i — ^4)] 



[A;i + 3A;4 — A;4A;6(3A;i + A;4)] 



^1^6 / , 2, 

g^2 v'*^4«^6 



ki) 



(66) 



(67) 



(68) 



The models represent an expanding, shearing but non-rotating universe 
in general. The models explode with a big bang at T = and the 
expansion in the models stops at T = 00. When k\ = Q then ^ = 0, 
which implies that = 0. Therefore, viscosity is due to expansion in the 
model. We take ki / 0. The spacetime is Petrov type I non-degenerate. 
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However, if k4 = 0, the spacetime reduces to Petrov type ID. For large 
values of T, the spacetime is conformally fiat. The expressions ^ and 
^ are found to be 



P 



kjkiikj - kl) + AAT^ 



(69) 
(70) 



The rate of expansion (Hubble parameters) in the direction of X, 
Y, Z are given by 

Hi = (71) 

kikeiki + k^) 



Ho 



H. 



2T 

kikQ{ki — k^) 

2r 



(72) 



(73) 



Since hmy^oo ^ 7^ Oj hence the models do not approach isotropy for 
large values of T. 



4. Solution in Absence of Shear Viscosity 

When 77 — 0, then the metric (27) leads to 

The pressure and density for the model (74) are given by 



1 



16M2T4 



A 



(M2 - L2) , 



(74) 

(75) 
(76) 



4.1. MODEL I: SOLUTION FOR ^ = 

When n = 0, Equation (31) reduces to ^ = ^0 (constant). Hence in this 
case Equation (75), with the use of (29) and (76), leads to 



87r(l + 7)p 



327rM2^or^ + - 



8M2T4 

Eliminating p{t) between (76) and (77), we obtain 
1 



(1 + 7)A 



i6M2r4 



QAttM^^oT"^ - (M2 - l2)(1 - 7) 



(77) 



(78) 
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4.2. MODEL II: SOLUTION FOR C = S.oP 

When n = I, Equation (31) reduces to ^ = ^op. Hence in this case 
Equation (75), with the use of (29) and (76), leads to 

— 

[2(l + 7)T^-eo]p=32^^ (79) 
EUminating p{t) between (76) and (79), we have 

[2(1 + 7)r2 - eo]A = [4r2 - 2(1 + 7)r2 + Co] (80) 

From Equations (78) and (80), we observe that the cosmological con- 
stant in both the models is a decreasing function of time and it ap- 
proaches a small positive value for large time (i. e. the present epoch), 
when — > 0, which is supported by the results from recent 
type la supernovae observations (Garnavich et al., 1998; Perlmutter et 
al.,1997, 1998, 1999; Riess et al, 1998; Schmidt et al, 1998). 
Some Physical Aspects of the Models: 

The scalar of expansion(^), shear (cr) and the non- vanishing components 
of conformal curvature tensor are given by 



2T2 



(81) 



2 _ (M2 + 3L2) 



48M2r4 



(83) 



The expressions | and ^ are as follows: 

^ = ^U(m2 + 3L2)5 (85) 
e 2V3M 

Prom Equation (85), it is observed that the shear is proportional to 
the expansion in the models. The models explode with a big bang 
at T = and the expansion in the models stops at T = oo. The 
spacetime is non-degenerate Petrov type I in absence of viscosity and 
if L = 0, then spacetime is Petrov type D. For large values of T, the 
spacetime is conformally flat. Since limT^oo f 7^ 0, hence the models 
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do not approach isotropy for large values of T. The metric (74) has 
singularity at T = 0. The models start from T = with big bang and 
goes on expanding till T = oo. 

5. Conclusions 

Some Bianchi type I anisotropic cosmological models with a viscous 
fluid as the source of matter are obtained. Generally, the models are 
expanding, shearing and non-rotating. It is observed that in the pres- 
ence of viscosity, the spacetime is Petrov type I D, whereas without 
viscosity the spacetime is conformally fiat. It is seen that the solutions 
obtained by Bali and Jain (1997) are particular cases of our solutions. 
In all these models, we observe that they do not approach isotropy 
for large values of time T. These models are new and different from 
those models obtained by Bali and Jain (1987, 1988) in which free 
gravitational field was assumed to be Petrov type D and non-degenerate 
for Marder (1958). 

The cosmological constant in all models given in sections 3.2 and 4 are 
decreasing function of time and they all approach a small positive value 
as time increases (i.e., the present epoch). The values of cosmological 
"constant" for these models are found to be small and positive which 
are supported by the results from recent supernovae la observations 
recently obtained by the High - Z Supernova Team and Supernova 
Cosmological Project ( Garnavich et al, 1998 ; Perlmutter et al, 1997, 
1998, 1999; Riess et al, 1998; Schmidt et al. , 1998). In section 3.1, 
our study is inconsistent with work done by Bali and Jain (1997). Our 
study for these models shows constant value of cosmological constant 
for large values of time and do not decrease with time i. e. universe is 
not expanding or may be steady state condition. Bali and Jain claim 
that the universe is expanding which does not match with our results. 
We are trying to find feasible interpretation to this case. Further study 
is in progress. 



Acknowledgements 

One of the authors (A. Pradhan) thanks to the Inter-University Cen- 
tre for Astronomy and Astrophysics, India for providing facility under 
Associateship Programmes where part of work was carried out. 



pumimalsty.tex; 7/02/2008; 3:00; p. 14 



SOME BIANCHI TYPE I VISCOUS FLUID COSMOLOGICAL MODELS 



15 



References 

Bali R. and Jain, V.C.: 1987, Astrophys. Space ScL, 139, 175. 
Bali R. and Jain, V.C.: 1988, Astrophys. Space ScL, 141, 207. 
Bali R. and Jain, V.C.: 1997, Astrophys. Space Sci., 254, 13. 
Bertolami, O.; 1986, Nuovo Cimento, B 93, 36. 
Bertolami, O.: 1986, FoHschr. Phys. 34, 829. 

Carroll, S.M., Press, W.H. and Turner, E.L.: 1992, Ann. Rev. Astron. Astrophys. 
30, 499. 

Chakrabarty, L, Pradhan, A. and Saste, N.N.: 2001, Int. J. Mod. Phys. D 10, 741. 
Dolgov, A.D.: 1983, in The Very Early Universe, eds. Gibbons, G.W., Hawking, 

S.W. and Siklos, S.T.C. (Cambridge University Press). 
Dolgov, A.D., Sazhin, M.V. and Zeldovich, Ya.B.: 1990, Basic Modem Cosmology, 

(Editions Frontiers). 
Dolgov, A.D. and Silk, J.: 1993, Phys. Rev. D 47, 3144. 
Dolgov, A.D.: 1993, Phys. Rev. D 48, 2499. 
Dolgov, A.D.: 1997, Phys. Rev. D 55, 5881. 

Garnavich, P.M. et al: 1998a, Astrophys. J. 493, L53, astro-ph/9710123; 

1998b, Astrophys. J. 509, 74, astro-ph/9806396. 
Heller, M. and Klimok, Z.: 1975, Astrophys. Space Sci. 33, 37. 
Johri, V.B. and Sudarshan, R.: 1988, Phys. Lett. A 132, 316. 

Kalyani, D. and Singh, G.P.: 1997, in New Direction in Relativity and Cosmology, 

eds. V. de Sabbata and T. Singh, Hadronic Press, USA, p. 41. 
Marder, L.: 1958, Proc. Roy. Soc. A 246, 139. 

Maartens, R.; 1995, Class Quantum Gravit. 12, 1455. 
Misner, C.W.: 1967, Nature 214, 40. 
Misner, C.W.: 1968, Astrophys. J. 151, 431. 
Murphy, G.L.: 1973, Phys. Rev. D 8, 4231. 
Nightingale, J. P.: 1973, Astrophys. J. 185, 105. 
Padmauabhan, T.: 2003, Phys. Rep. 380, 235; hep-th/0212290. 
Padmanabhan, T. and Chitre, S.M.: 1987, Phys. Lett. A 120, 433. 
Pavon, D., Bafaluy, J. and Jou, D.: 1991, Class. Quant. Grav. 8, 357; 

1996, "Proc. Hanno Rund Conf. on Relativity and Thermodynamics", Ed. S. D. 

Maharaj, University of Natal, Durban, p. 21. 
Perlmuttor, S. et al: 1997, Astrophys. J. 483, 565, astro-ph/9608192; 

1998, Nature 391, 51, astro-ph/9712212; 

1999, Astrophys. J. 517, 565, astro-ph/9608192. 

Peebles, P.J.E. and Ratra, B.: 2003, Rev. Mod. Phys. 75, 559, astro-ph/0207347. 

Pradhan, A., Sarayakar, R.V. and Bccsham, A.: 1997, Astr. Lett. Commun. 35, 283. 

Pradhan, A., Yadav, V.K. and Chakrabarty, I.: 2001, Int. J. Mod. Phys. D 10, 339. 

Pradhan, A., V.K. Yadav and Saste, N.N.: 2002, Int. J. Mod. Phys. D 11, 857. 

Pradhan, A. and Yadav, V.K.: 2002, Int. J. Mod. Phys. D 11, 839. 

Pradhan, A. and Aotcmshi, I.: 2002, Int. J. Mod. Phys. D 11, 1419. 

Pradhan, A. and Pandcy, H.R.; 2003, Int. J. Mod. Phys. D 12, 941. 

Pradhan, A. and Pandey, O.P.: 2003, Int. J. Mod. Phys. D 12, 1299. 

Pradhan, A. and Singh, S.K.: 2004, Int. J. Mod. Phys. D 13, 503. 

Pradhan, A. and Srivastav, S.K. and Jotania, K.R.: 2QQA,Czech. J. Phys. 54, 255. 

Pradhan, A. and Yadav, L. and Yadav, A.K.: 2004, Czec/i. J. Phys. 54, 487. 

Ratra, B. and Peebles, P.J.E. : 1988, Phys. Rev. D 37, 3406. 

Riess, A.G. et al: 1998, A.stron. J. 116, 1009, astro-ph/9805201. 

Roy, S.R. and Prakash, S.: 1977, Ind. J. pure Appl Math. 8, 723. 



pumimalsty.tex; 7/02/2008; 3:00; p. 15 



16 



A. PRADHAN AND P. PANDEY 



Sahni, V. and Starobinsky, A.: 2000, Int. J. Mod. Phys. D 9, 373, gr-qc/9904398. 
Santos, N.O., Dias, R.S. and Bancrjoo, A.: 1985, J. Math. Phys. 26, 878. 
Schmidt, B. P. et al: 1998, Astrophys. J. 507, 46, astro-ph/9805200. 
Singh, T., Beesham, A. and Mbokazi, W. S.: 1998, Gen. Rel. Grav. 30, 537. 
Tsagas, C.G. and Maartens, R.: 2000, Class. Quant. Grav. 17, 2215, gr-qc/9912044. 
Vishwakarma, R. G. and Abdussattar: 1999, Phys. Rev. D 60, 063507. 
Vishwakarma, R.G.: 2000, Class. Quant. Crav. 17, 3833. 
Vishwakarma, R.G.: 2001, Class. Quant. Crav. 18, 1159. 
Vishwakarma, R.G.: 2001, Gen. Rel. Grav. 33, 1973. 
Vishwakarma, R.G.: 2002, Mon. Not. R. Astron. Soc. 331, 776. 
Vishwakarma, R.G.: 2002, Class. Quant. Crav. 19, 4747. 
Weinberg, S.: 1971, Astrophys. J. 168, 175. 

Weinberg, S.: 1972, Gravitation and Cosmology, Wiley, New York. 
Zeldovich, Ya. B.: 1968, Sov. Phys.-Uspekhi 11, 381. 
Zimdahl, W.: 1996, Phys. Rev. D 53, 5483. 



purnimalsty.tex; 7/02/2008; 3:00; p. 16 



